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In a recent paper it was proposed that for some nonlinear shell models of turbulence one can 
construct a linear advection model for an auxiliary field such that the scaling exponents of all 
the structure functions of the linear and nonlinear fields coincide. The argument depended on an 
assumption of continuity of the solutions as a function of a parameter. The aim of this paper is to 
provide a rigorous proof for the validity of the assumption. In addition we clarify here when the 
swap of a nonlinear model by a linear one will not work, 



I. INTRODUCTION 

Shell models of turbulence [l], Q serve a useful purpose 
in studying the statistical properties of turbulent fields 
due to their relative ease of simulation. In particular, 
shell models allowed accurate direct numerical calcula- 
tion of the scaling exponents of their associated structure 
functions, includingconvincing evidence for their univer- 
sality 0,a,!,!,0,ii- 

In contrast, for the Navier-Stokes 
equations that model actual fluid turbulence simulations 
are very much harder, and in addition one still does not 
know whether these equations in 3-dimcnsions are math- 
ematically globally well posed. This problem does not 
exist in shell models f2lj . adding to their numerical at- 
tractiveness a possibility to prove various properties and 
results rigorously [2l|, [22] ■ 

Consider for example the Sabra shell model Q which, 
like other shell models of turbulence, is a truncated de- 
scription of the dynamics of Fourier modes, preserv- 
ing some of the structure and conservation laws of the 
Navier-Stokes equations: 



(-r + vk 2 n )u n = i[k n+ iu* n+1 u n+2 - 5k n u* n _ 1 u n+ i 
+ (l-S)k 

n — 1 ^ n — 1 Uji — 2 } + In ■ (1) 

Here u n , with n = 0, 1, 2, . . . and the boundary condi- 
tions w_2 = U-i = 0, are the velocity modes restricted 
to 'wavevectors' k n = kofi n with fco determined by the in- 
verse outer scale of turbulence. The model contains one 
additional parameter, S, and it conserves two quadratic 
invariants (when the force and the dissipation term are 
absent) for all values of 5. The first is the total en- 
ergy ~^2 n \u n \ 2 and the second is ^ n {~^-) n k c ^\u n \ 2 , where 
a = log M (l - 5). 

The scaling exponents are properties of the structure 
functions. To define the structure function we introduce 
an average over time according to 



In practice, however, we take 
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(A(t)) = A(t)dt , 



(3) 



(A(t)) 



1 f T 
km - / A(t)dt 



(2) 



for some T large enough, but finite. Our rigorous re- 
sults also refer to this definition of an average over time. 
For values of the viscosity v small enough, and for a suf- 
ficiently large amplitude of the random force /„ there 
exists a large range of values k n where numerical simu- 
lations show that structure functions follow a power-law 
behavior. The low-order structure functions and the as- 
sociated scaling exponents are 

S 2 (k n ) = {u n u* n ) ~ k~ C2 , (4) 
S 3 (k n ) = 3(u„_iu„u* +1 ) ~ k~ C3 , (5) 
etc. for higher order S p {k n ) ~ k~^ p . 

The values of the scaling exponents were determined 
accurately by direct numerical simulations. Besides £3 
which is exactly unity @ , all the other exponents £ p ap- 
pear anomalous, differing from p/3. Numerical evidence 
is that the scaling exponents are also universal, i.e. they 
are independent of the forcing /„ as long as the latter is 
restricted to small n • It was shown that for < S < 1 
the leading scaling exponents are determined by the cas- 
cade of the energy invariant from large to small scales. 
For 1 < S < 2 the second invariant \u n \ 2 be- 

comes important in determining the leading scaling ex- 
ponents of the structure functions of u n . In the bulk of 
this paper we consider the situation with < 5 < 1, but 
return to the other case in Sect. IIV1 in order to clarify the 
role of invariants in determining the scaling properties. 

In a recent publication further insight to the anomaly 
of the exponents of the nonlinear problem (for the field 
u n ) was sought by relating them to the scaling exponents 
of a linear model for a field w n . The linear model was 
constructed such that its scaling exponents should be the 
same as those of the nonlinear problem. The equations 
for this field are constructed under the following require- 
ments: (i) the structure of the equations is obtained by 
linearizing the nonlinear problem and retaining only such 
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terms that conserve the energy; (ii) the resulting equa- 
tion is identical with the sabra model when w n = u n ; (hi) 
the energy is the only quadratic invariant for the passive 
field in the absence of forcing and dissipation. These 
requirements lead to the following linear model: 

dw^ i , . o , , 

= W) - Vk n W n + f n , (6) 

where the advection term is defined as 

w) = fc n+ i[(l + 8)u n+ 2W* n+1 + (2 - &)u* n+1 w n +2\ 
+k n [(l - 28)u* n _ 1 W n +l - (1 + S^n+lW*^] 
+fc„_l[(2 - S)u n -iW n -2 + (1 - 2S)u n - 2 W n -l] , (7) 

together with u_i = u_2 = W-i = w_ 2 = 0. Observe 
that when w n = u n this model reproduces the Sabra 
model, and also that the total energy is conserved be- 
cause J2 n ^[^n(u, w)w^] = 0. The second quadratic 
invariant is not conserved by the linear model (J6j) . Fi- 
nally, both models have the same 'phase symmetry' in the 
sense that the phase transformations u n — > u n exp (i(j) n ) 
and w n — y w n exp (i6 n ) leave the equations invariant iff 
4> n -i + 4> n = 4> n +i, n -i + &n = n +i- This identical 
phase relationship guarantees that the non- vanishing cor- 
relation functions of both models have precisely the same 
forms. Thus for example the only second and third cor- 
relation functions in both models are those written ex- 
plicitly in Eqs. g]) and ©. 

The advantage of the linear model is that correlation 
functions are advanced in time by a linear propagator 
0, EH, Ell- The linear model possesses "Statistically 
Preserved Structures" (SPS) which are evident in the 
decaying problem Eq. with /„ — 0. These are left 
cigenfunctions of eigenvalue 1 of the linear pr opag ators 
for each order (decaying) correlation function [14j . For 
example for the second order correlation function denote 

( 2^ 

the propagator P^ ' n , (t|io); this operator propagates any 
initial condition (w n w^)(to) (with average over initial 
conditions, independent of the realizations of the advert- 
ing field u n ) to the decaying correlation function (with 
average over realizations of the advecting field u n ) 

= ^(*l*o)K'0(*o) • (8) 

(2) 

The second order SPS, Z„ , is the left eigenfunction with 
eigenvalue 1, 

Z^=Z^Pi%(t\t ). (9) 

(2) 

Note that Z n is time independent even though the oper- 
ator ' n , (t\ta) is time dependent. Each order correlation 

function is associated with another propagator (t\to) 
and each of those has an SPS, i.e. a left eigenfunction 
of eigenvalue 1. These non-decaying eigenfunctions 
scale with k n , ~ k n , and the values of the ex- 

ponents £ p arc anomalous. Finally, it was shown that 
these SPS are also the leading scaling contributions to 




iog(k n ) 

FIG. 1: The sixth order structure function of the field w n 
in Eqs. CEO} for A = 10 -1 , 1(T 3 and 1CT 5 , together with the 
sixth order structure function for the Sabra model <jTJ and 
for the linear model (|6]), respectively. The structure function 
of the field u n for A > are not shown since they are indis- 
tinguishable from those of the w n . Inset: log-log plot of the 
fourth-order correlation function F2, 2(^71,^7) vs. k n calcu- 
lated for the linear field (+) and for the nonlinear field (solid 
line) at A = 0. 

the structure functions of the forced problem ([6]) [3, [3 ■ 
Thus the scaling exponents of the linear problem 
are independent of the forcing f n , since they are 
determined by the SPS of the decaying problem. 

To connect the linear model to the nonlinear problem 
one considers the system of two coupled equations 

du n i iX 2 

~dT = 3*"( M ' u > + y ®n{w, u) - vk n u n + f n , 

^ = i$ n («, + y $„(«;, to) -i/fc> n + / n (10) 

with A being a real parameter and f n and /„ being dif- 
ferent realizations of the same random force. Observe 
that for any A ^ 0, the two equations in (fTU|) exchange 
roles under the change Xw n <-> u n . Thus if the scaling 
exponents £ p and Cp of the two fields exist (i.e. a true 
scaling range exists), they must be the same for all A 7^ 0. 
For A = we recover the equations for the nonlinear and 
the linear models, Eqs. |T]) and ([5]). In [|o| it was as- 
sumed that the scaling exponents of either field exhibits 
no jump in the limit A — > 0. Indeed, Ref. [2(| also pre- 
sented strong evidence for the validity of this assumption 
(see Fig. [TJ), but no mathematical proof was provided. 

The aim of this paper is to close this gap. The main 
result of our paper states that the solutions of the system 
pop exist globally in time and depend continuously on 
the parameter A, including the limit A — > 0. In particular 
we will show that if the structure function of u„ and w n 
exhibit the same scaling exponents for any A ^ 0, they 
will also have the same scaling exponents in the limit 
A — » (including A = 0). We would like to stress here 
that our rigorous results arc valid for the structure func- 
tions, calculated over large, but finite, fixed interval of 
time, which is consistent with the definition ([3]) of the 
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long time average. In addition, we would like to note 
that the numerical results correspond to the equations 
with the stochastic implementation of the forcing, while 
our rigorous proofs deal with the deterministic force that 
depend on time. The statement and the proof of the 
main theorem will be given in Section IIIII The proof is 
based on the results on the global existence and unique- 
ness of solutions of Eq. ((T|) , obtained previously in [21j . 
In the following section we present the necessary mathe- 
matical definitions and formulate the essential statement 
from [2l|. 

II. FUNCTIONAL SETTING AND PREVIOUS 
ANALYTIC RESULTS 

A. Functional Setting 

Following the classical treatment of the Navier-Stokcs 
and Euler equations we re- write the sabra shell model for 
an infinite vector u = (uq, u\, . . .) 

du i T , . , , 

— + vAu + -*(«,«) = / , (11) 

together with the initial conditions u(t = 0). Intro- 
duce a Hilbert space H which is the space of infinite 
vectors equipped with the scalar product (•,•) and the 
corresponding norm | • | defined as 

co / oo s 1/2 

(u,«) = 5^u n < J M = (]TK| 2 j , (12) 

n=0 ^ n=0 ' 

for every u,v £ H. The space H is a space of all the 
velocity vectors having finite energy. 

The linear operator A, with a domain D(A) dense 
in H, is a positive definite, diagonal operator defined 
through its action on the elements u by 

Au = (kino, kfui, . . .) (13) 

where the eigenvalues fc 2 satisfy k n = k /j, n . Using the 
fact that A is a positive definite operator, we can define 
the powers A s of A for every s £ TZ 

Vit = (u , ui, u 2 , ■ • ■), A s u = (/c 2s m , k 2s ui, kj s u 2l ■ ■ ■)■ 

(14) 

The space D{A S / 2 ) is the domain of the operator A s l 2 
and we denote 

oo 

V s = D{A S ' 2 ) = {u = (uo.ui, ua, ■ ■ ■) *i*Kf < 

3=0 

(15) 

which arc Hilbert spaces equipped with the scalar prod- 
uct 

(u, v) s = {A s ' 2 u, A s ' 2 v), Vu, v £ D{A S / 2 ), (16) 

and the norm \u\ 2 = (u,u) s , for every u £ D(A S ^ 2 ). 
Since V s contains velocity vectors with s "derivatives" , 

V a C V = H C V- s , Vs > 0. (17) 



The case of s = 1 is of a special interest for us. We 
denote V = Z^A 1 / 2 ) a Hilbert space equipped with a 
scalar product and a corresponding norm 

((«,«)) = {A 1 ' 2 u 1 A 1 ' 2 v), \\u\\ 2 = ((«,«)), (18) 

for every u, v £ V. In addition the following relation 
holds 

fciM < ll«l| ■ (19) 

In what follows we will need the interpolation inequality 
for the spaces V s . 

Lemma 1. Let s > 0, then for all u £ V s and < 

s' < s 

\u\ s , < Itif-'V-M;'/- . (20) 

Proof. The Lemma follows by a simple application of 
Holder inequality. ■ 

The bilinear operator |*(tt, w) is defined above, cf. 
Eq. together with # = ($ , $i, . . .). In [2l| it 
was shown that such a definition of the bilinear oper- 
ator makes $ and element of H whenever u £ H and 
w £ V. For any u,v £ H and w £ V one proves [2l| 

| (*(«,«), to) | < C\u\ \v\ \\w\\ . (21) 

for some positive constant C. In addition the conserva- 
tion law is written is the present notation as 

3f(*(u,io),io) = . (22) 

All our rigorous results are valid for the deterministic 
forcing /. Therefore, in order to account for the stochas- 
tic implementation of the forcing term we will assume 
that / depends on time, but always stays bounded in the 
if- norm. More precisely, define the space L°°([0,T), H), 
for some < T < oo, as a space of functions of the time 
interval [0, T) with values in H and the norm defined as 

H/lloo = sup \f(t)\ . (23) 

0<t<T 

In what follows, we will assume that the forcing term / 
satisfies / £ L°°([0,oo),H). 

B. Summary of Previous Results 

In [2l| Eq. (TTJ) was studied, and the relevant results 
can be formulated as the following theorem: 

Theorem (CLT 06): The solution of Eq. ((J) exists 
globally in time for any initial condition u(0) in H, and 
for any / in L°°([0, oo), H). Moreover, the solutions are 
unique and the energy of the solution u(t) is bounded for 
all times: 

sup \u(t)\ 2 < K (k ,(i,v,5,u(0),f) , (24) 

0<t<T 
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where the a-priori constant Kq depends on all the param- 
eters of the equations, on the forcing and on the initial 
conditions. 

If in addition we assume that the forcing / = 
(/o,/i,...) acts on the finite number of shells, namely, 
if there exists N > 0, such that /„ = 0, for all n > N, 
then the solution u = (uq, Ui, . . .) has an exponentially 
decaying spectrum | function of k n , and in par- 

ticular 

sup |«(t) |f < «T B (*o, ^ i/, <5, «(0), /), (25) 

vkl<t<T 

for any s > and the a-priori constants K s depend on 
all the parameters of the equations, on the forcing and 
on the L°°([0, oo), ff)-norm of the initial conditions (see 
definition J53J). 



III. THE MAIN RESULT 

The main statement of this paper is that the system 
of equations (JTOJ) is globally well posed for all real A, and 
that the solutions depend continuously on the parameter 
A. In particular, as A — > the solution of the system 
converge uniformly on any finite interval of time to the 
solution of the system with A = 0. This statement is 
formulated as follows: 

Theorem 1. Let it(0; A), w(0; A) be given in H, 
and the forces /,/ in L°°([0, T],H). Denote by 
u(t; X),w(t; X) the solutions of the coupled system (fTQ|) . 
Then the following holds 

1. The solution of the system ([Hi]) exists globally in 
time for every A. Moreover, energy of the solutions 
|it(£;A)| 2 and |«>(i;A)| 2 are bounded for all times, 
where the bounds depend on all the parameters of 
the equations, A, the forcing and the initial condi- 
tions. 

2. For any T > 0, 

lim sup \u(t; A) — u(t; 0)| 2 = 0, 

*->°0<t<T 



and 



lim sup \w(t; A) - w(t; 0)| 2 = 0, 
A - >0 o<t<T 



where u(t;0),w(t;0) are the solutions of (fTU|) with 
A = 0. 

3. If in addition we assume that the forces /, / act 
only on the finite number of shells then for any 
T > vk 2 , and s > 0, 

lim sup \u(t; A) - u(t; 0)| 2 = 0, 

vkl<t<T 



and 



lim sup \w(t; A) - w(t; 0)| 2 = 0, 



^— >0 ,,1.2 



uki<t<T 



where u(t; 0), w(t; 0) are the solutions of (fTU|) with 
A = 0. 



Proof. Part 1. of this theorem follows from defining a 
new variable q(£; A) = u(t; A) + Xw(t; A). This variable 
satisfies Eq. ([I]) with a forcing f + Xf. Accordingly theo- 
rem CLT06 provides the proof that q(t; A) exists globally 
in time for every A. Next observe that the system of 
equations (|10|) can be rewritten in the form 



~dt 

dw„ 

~~dT 



$ n (q,u)) - vk 2 n u n + f n , 

w) - vk 2 n W n + f n . 



(26) 



This form of writing shows that the the fields u and w 
satisfy linear diffusion advection equations advected by 
a smooth field q. Accordingly both fields remain smooth 
for all time and all A. In addition, using the relation (|2"4"|) 
one is able to derive the bounds for the energies of the 
solutions. 

To prove Part 2. of the theorem we need first 
Proposition 1: Denote by q(t;0) the solution of the 
Sabra model (JTJ) with initial data q(0; 0) in H. That is 
also the solution of the first of equations (fTDj) when A = 0. 
Then 



lim sup \q(t; A) - q(t;0)\ = 0. 

0<KT 



Moreover, if /, and / act only on finite number of shells, 
then for any s > we have 

lim sup | g (t;A)-g(t;0)| a = 0. 



vkj,<t<T 



Proof. Let us denote by Aq = q(t; A) — q(t; 0). Then 
Aq satisfies the equation 

^ + vAAq - l -<t>(q(t; A), Aq) - ^*(Aq, q(t; A)) 



+-<f>(Aq,Aq) = Xf 
Aq(0) = Xw(0). 



(27) 
(28) 



Take now the inner product in H of the above equation 
with Aq. Computing the real part and using Eq. (|2"2"1) 
we find 

~| Aq\ 2 +u\ | Aq\ | 2 -5R(^*(Aq, q(t; A), Aq) = R(A/, Aq) 

(29) 

Using the Cauchy-Schwarz inequality and the relation 
(12T|) we get 



ld_ 

2dt 



Aq| 2 + H|Aq|| 2 < 



< |(-*(Aq,q(t;A),Aq)| + |(A/,Aq)| 

<C|Aq| ||Aq|| |q(i;A)|+A|/| |Aq| , (30) 

Applying the Young's inequality (ab < a 2 /2 + & 2 /2) twice 
and using the inequality (|19p . we have 
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l -±\Aq\ 2 + v\\Aq\\ 2 <^-\Aq\ 2 



*fcA)|' + ^|/ 



2 + |HAq|| 2 . 



r 



Using the fact that the relation holds for any t > 
and the fact that / e L°°([0, oo),H) (see definition 1(25)1). 
we obtain 



d_ 

dt 



9 ft^r 12 

i Ag | 2 <^-|A g | 2 + ^n/r 2 



2A 2 
i/fc 2 ' 



By Gronwall's inequality we conclude 
|A 9 (i)| 2 < 

2K 2 C 2 A 2 llfll 2 2K n° 2 

< eH^*| A q(0)| 2 + y^ffi (e -B-« - 1) 



A 2 C 2 fc 2 

1 A 2 llfll 2 2K n° 2 

-*A 2 M0)| 2 + HJll~^-S-* 



A 2 C 2 fc 2 



1). (31) 



Therefore, for any T > 0, 



lim sup |Aq(t)| 2 = 0, 

A ^°0<t<T 

and the first statement of the proposition follows. 



The second statement of the proposition follows from 
the boundedness of q(t; A) and q(t; 0) in higher order 
norms after a short transient period of times, say fk^, 
(provided that the forces /, and / act on the finite num- 
ber of shells as required by Theorem CLT06) and the 
interpolation inequality ((20)) . ■ 

Finally, we are ready to finish the proof of the main 
theorem. Let us fix s > 0, and show that 

lim sup \w(t; A) - w(t; 0)| 2 = 0. (32) 

ukl<t<T 

Denote Aw — w(t; A) — w(t; 0). The function Aw satis- 
fies the equation 

j t Aw + vAAw - l -$(q{t- A), w(t; A)) 



+ -*(g(t;0),t»(i;0)) = 0. 



(33) 



We rewrite it in the form 



—Aw + vAAw - -*(<?(*; A) - q(t; 0), w(t; 0)) - -*(qr(*; A), Aw) = 0, 



(34) 



and as before, taking the inner product in H with Aw, computing the real part and using Eq. (|22p . we obtain 

1 d 



lt t Aw\ 2 + v\\Aw\\ 2 - »(-*(g(t; A) - q(t; 0), w(t; 0)) , Aw) = 

Applying, subsequently the Young's inequality and the inequality (|21j) . we get 
1 d 



2 dt 



A™| 2 + H|AH| 2 < |(-#(g(f;A)-g(t;0),ti;(t;0)),At»)| < 



< C\\Aw\\ \w(t; 0)| \q(t; A) - q(t; 0)| < ^-\w(t; 0)| 2 \q(t; A) - q(t; 0)| 2 + ^\\Aw\\ 2 . 



r 



(35) 



of the main theorem. It follows from Proposition 1 that 
sup 0<t<T \q(t; A) — q(t;0)\ 2 — > 0, as A — > 0. Hence we 
may conclude that 

lim sup \Aw(t)\ 2 = 0. 

0<t<T 

By virtue of the interpolation inequality ((20)) one can fol- 
low steps as in the proof of Proposition 1 to show that the 
where we used the fact that Aw(0) = and |u;(i;0)| 2 < higher order norms of Aw also converge to uniformly 
M, for some constant M, as was stated in the Part 1 in the time interval vk'k < t < T, as A — > 0. To finish the 



It follows that 

±\Aw\ 2 < — \w(t;0)\ 2 \q(t;\)-q(t;0)f 
dt v 

and integrating over (0, t) we conclude 



\Aw{t)\ 2 < 



C 2 M 



t\q(t;\)-q(t;0)f 



shell numeber n 



shell numeber n 



FIG. 2: Double logarithmic plot of the second order structure 
functions (u„u*) and {w n w^} of the system of equations (|10[) 
for S = 1.25 and A = 0. 
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FIG. 3: Double logarithmic plot of the second order structure 
functions («„«*) and (w n w^} of the system of equations (fT0|) 
for S = 1.25 and A = 1. 

proof, on can observe that 

\u(t; A) - u(t; 0)| = |g(i; A) - q(t; 0) - Aw(t; A)| 
<\q(t;\)-q(t;0)\ + \\w(t;X)\, (36) 

and from the fact that |to(t;A)| is globally bounded in 
time it follows from Proposition 1 that 

lim sup \u(t;X)—u(t;0)\=0. 

A ^0 <t<T 



A. Consequences for Structure Functions 

Using the form of Eqs. (f2"6"| we conclude that for f n and 
/„ being different realizations of the same random force, 
whenever the structure functions of u n and w n exhibit 
the same scaling exponents for all finite A, they must 
also exhibit the same scaling exponents for A = 0. 




5 10 15 20 25 

shell numeber n 



FIG. 4: Double logarithmic plot of the second order structure 
functions {u n Un) and {w n Wn) of the system of equations (|10p 
for S = 1.25 and A = 0.1 (upper panel) and A = 0.01 (lower 
panel). 



IV. WHEN CAN THE NONLINEAR MODEL 
EXHIBIT SCALING EXPONENTS THAT ARE 
DIFFERENT FROM THE LINEAR MODEL 

In this section we turn our attention to situations when 
the nonlinear and the linear models cannot exhibit the 
same scaling exponents. The theorem as stated and 
proven still holds, but as we shall see this is a situa- 
tion for which the two fields u n and w n cannot exhibit 
the same scaling exponents for all A ^ 0. A case in 
point is the same nonlinear Sabra model with 1 < 5 < 2. 
In Fig. [2] we show the second order structure functions 
(w n u*) and (w n w^) obtained from simulating Eqs. (jTUJ) 
for A = and S = 1.25. As expected, the scaling of w n 
is influenced by the cascade of energy, whereas that of 
u n by the cascade of the second invariant. As a result 
the scaling exponents are distinctly different. The same 
system of equations for A = 1 is symmetric in w n and 
u n . Indeed, in Fig. [3] we show the result of simulations 
for A = 1, for which the second order structure function 
of the two fields is identical. Now however we cannot 
expect that this identity will persist for A — > 0. In Fig, 
13] we show the results of simulations for the same system 
of equations for A = 0.1 and A = 0.01. 

To understand the results of the simulations we note 
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that when A ^ the second invariant of the equation for 
u n is destroyed, and one could think that the scaling of 
u n should be dominated by the energy invariant. This is 
certainly true for A = 1. But now when A decreases to- 
wards zero, we should reconsider the system of equations 
(|10|) by renaming w n = Xw n ■ Substituting this re- named 
variables into Eqs. (|10[) and re-arranging, we read 

—rr = ^n(u,u) + ?-<S> n (w,u) - vk*u n + f„ , 
at 3 3 

~df = 3* n ( u '*) + ^n(w,w) ~ "K^n + A/ n (.37) 

Thus the net result of the transformation is that the equa- 
tions for u n and w n are the same, but the forcing of w n 
becomes weaker and weaker as A — » 0. Accordingly, while 



the second invariant is still destroyed as a true invariant 
for any value of A, for small A the strength of the term 
|$„(w,u) diminishes, allowing a cross-over behavior in 
the scaling of u n , precisely as we see in Fig. @] 
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